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Abstract. We propose the use of the Dirac-Born-Infeld action in the phenomenological description of 
graphene sheet dynamics and interactions. Both the electronic properties of the Dirac fermions and the 
overall dynamics can be incorporated into this model. Classical static configurations, as well as quantum 
fluctuations of the membrane degrees of freedom can be studied in this framework. This makes it an 
interesting tool for Casimir physics and novel QED processes. 



1 Introduction 

The Dirac-Born-Infeld (DBI) action [fll], originally writ- 
ten down as a way of extending classical electrodynam- 
ics and addressing the electron's self energy, is an essen- 
tial tool if one wants to consider the dynamics of higher 
dimensional objects e.g. strings and membranes. Its use 
in string theory has been especially important for the 
case of Dirichlet-branes (D-branes) [3j|4]. These are non- 
perturbative extended objects which have the ends of fun- 
damental open strings ending on them. At the same time, 
the low energy dynamics of the brane are described sim- 
ply by a gauge theory living on this surface. There are 
two interesting features we wish to draw attention to here. 
Firstly, the degrees of freedom of the brane couple natu- 
rally to the fields in the bulk. Secondly, the microscopic 
theory can be used to calculate perturbatively the cou- 
pling constants between the brane and bulk fields. A nat- 
ural question arises as to whether this type of dynamics 
can be applied to other systems. 

Graphene has the interesting property that the degrees 
of freedom describing the low energy electronic proper- 
ties are relativistic Dirac fermions [5,6,7,8 . This together 
with its one-atom thickness allow in principal for novel ap- 
plications. The Dirac action allows one to calculate phys- 
ical observables such as the electrical conductivity or the 
density of states. If we now want to consider the larger 
dynamics of the graphene sheet, not just the electronic 
properties, we need to be able to describe the graphene 
membrane as a whole [9] . Macroscopic graphene sheets are 
available pUUTT] that are in principle dynamical objects in 
their own right. This together with the tension and tor- 
sion p!2l[T3] necessitates the need for a suitable dynamical 
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theory that captures the relevant field theory aspects in a 
unified fashion. 

In this paper we consider using the DBI action (and 
suitable generalisations thereof) as a phenomenological 
tool for studying the low energy dynamics of graphene 
sheets. In particular, the low energy sector should not be 
able to resolve the hexagonal lattice structure of graphene 
(for perturbative field theory calculations the lattice struc- 
ture will act as an ultraviolet cutoff), whereby collective 
coordinates can simply be written down. One novel and 
new feature for this situation is the presence of real fermionic 
degrees of freedom - the quasi Dirac particles require the 
DBI action to be modified for their inclusion. This may 
be seen as a completion of the DBI action (which usually 
involves only bosonic fields) that has physical relevance 
and can be performed in a manner outlined in |14) . The 
low energy description of the hexagonal lattice (i.e. its 
position in space) trades the lattice structure and the in- 
dividual nuclei simply for a world sheet coordinate and 
an embedding. In a similar fashion, the quasi Dirac parti- 
cles become the fermionic completion of this embedding. 
One final factor that can be included is the coupling to 
the bulk electromagnetic field. This is incorporated as a 
polarisation type term that lives on the membrane itself, 
as well as from a covariant derivative of the fermions. We 
shall then work out some of the consequences due to the 
nature of their mutual couplings. 

The outline of this paper is as follows. In Section [5] 
we describe the basics of the DBI action to fix notations 
and conventions. In Section [3] we generalise this appro- 
priately to include the Dirac fermions on the membrane 
and possible couplings between the membrane itself and 
the bulk electromagnetic field. In Section |4] we perform an 
evaluation of the DBI action for a dielectric background 
that has only a one dimensional dependence and a sim- 
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pie graphene drum placed in a constant magnetic field. 
In Section [5] an expansion of the fermionic part of the 
action is performed to obtain a modified Dirac equation 
due to the presence of the polarisation term. This paves 
the way for three calculations - the spectrum of Landau 
levels, an induced Chern-Simons term and finally a modi- 
fied Schwinger effect. In Section [6] a path integral picture 
is presented whereby a static potential between between 
two parallel graphene membranes is obtained by consid- 
ering the S- matrix. Finally in Section [7] we draw our con- 
clusions. 



2 Basics of the DBI action 

To fix our notation and conventions, we briefly describe 
what the DBI action consists of in the relativistic setting 
of fleld theory. Consider ad — l-f2 membrane (two spatial 
dimensions, one time) embedded in the ambient D = 1-1-3 
bulk spacetime (three spatial dimensions, one time, see 
Figure [T|). The DBI action is then essentially the area of 
this space (denoted by the manifold ^73). One such model 
takes the form 



S3^n d3^V|det[3ab+J-ab]|. 



(1) 



S3 



The quantities and conventions used here are: T3 is the 
tension (or mass per unit area) of the membrane; ^'^ are 
the curved world sheet coordinates of the membrane with 
index a = 0,1,2; the embedding coordinates of the mem- 
brane are X^'^ {£,"') in the bulk space with the index AI — 
0, 1,2,3; and a bulk metric Gmn{X) field [4 . With this 
structure the bulk metric is said to be pulled back to 
the world sheet with the relation gab = GMNdaX^dtX^ 
and we can now calculate the area of the membrane. We 
have also included a d-dimensional U{1) gauge field on 
the membrane given by Aa, leading to the field strength 
J^ab = idA)ab (which is invariant under a gauge transfor- 
mation). This is a simple illustration of how one can also 
include into the action fields that live on the membrane. 
We will do this in particular for the Dirac fermions that 
are to describe the electronic properties of the membrane 
as well as for the electromagnetic field strength itself after 
performing an appropriate pull-back. 

Using the reparametrisation invariance of the DBI ac- 
tion, one can choose a gauge (sometimes called the static 
gauge or the Monge gauge) where the world-sheet coordi- 
nates coincide with the physical coordinates of the mem- 
brane leaving one coordinate orthogonal to the membrane. 
This describes the embedding of the membrane in the bulk 
space:- 



x^iO = C^ x^ ^ zix°,x\x^). 



(2) 



With this choice of gauge one sees that the DBI action 
takes the form of a non-linear field theory on the mem- 
brane. There is one scalar field that describes the embed- 
ding of the membrane in the bulk space. One can consider 



Local world-sheet coordinates 



Graphene membrane 
world-sheet Ej 




Generic fields 
defined oni^j 



Fig. 1. The d = 1 + 2 membrane embedded in the D = I + 
3 bulk space. The coordinates ^^ are the (in general) curved 
world sheet indices of the membrane, whilst the X^^ are the 
curved bulk coordinates. The membrane is then given as the 
embedding X'^' {^). Generic fields 'PatiO can be defined on the 
membrane. 



fluctuations around this embedding by expanding the de- 
terminant and the square root for a particular background 
metric. The expansion then gives a fully interacting field 
theory as a derivative expansion which can then be quan- 
tised. It also affords a description of minimum energy sur- 
faces for static configurations of membranes that are sub- 
jected to some potential via the bulk metric. 



3 The graphene DBI action 

We now seek to write down a phenomenological version 
of the DBI action for a graphene membrane. The degrees 
of freedom of the membrane are coupled to a set of bulk 
fields providing the necessary interactions. This is natu- 
rally a low energy description where the hexagonal lattice 
structure of graphene is simply replaced with the embed- 
ding coordinates described above, just as the fermions are 
replaced with the Dirac spinor fields. The underlying lat- 
tice structure also provides an UV cutoff for any of the 
effective field theory correlation functions that are calcu- 
lated on the membrane. 

Indeed in [15], the field equation for the bulk photon 
is written down with the addition of a static defect term 
to account for the interaction with the graphene fermions; 
it takes the form 



OmF^''' +6{X^-zo)n 



MN 



A 



M 



0. 



(3) 



Here 77*^^ is the one-loop polarisation operator for the 
photon on the graphene sheet extended to D = (1 + 3) 
dimensions and zq is the position of the graphene sheet 
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in the X^ direction. As it stands the graphene sheet is a 
static object that interacts only via its fermions. Our task 
is now to make the whole object dynamical by construct- 
ing a suitable low energy DBI action. 

The degrees of freedom on the membrane are the mem- 
brane coordinates X^{£^°') themselves, together with the 
d = 1 + 2 spinors '0('C") that describe the fermion-hole 
pairs on the graphene sheet itself. Both of these interact 
with the bulk 13 = 1 + 3 electromagnetic potential Am- 
The first quantity to form is a spinorial covariant deriva- 
tive on the membrane 

Dai' = a,V - eAMdaX^'i: - ]u:^^[n, r,]ij, (4) 

where Ua is a one-form spin connection (compatible with 
the pull-back metric gab), and F"" are the d dimension 
gamma matrices that are rescaled with the Fermi velocity 
vf as in |15| . The inclusion of torsion terms as discussed 
in [T^ is also possible should it need to be present. The 
membrane phonons, the bulk photon field and the curva- 
ture of the membrane can be incorporated in this way. In 
graphene there are four species of fermions, so we should 
really include an extra index on the spinor to account 
for this; however we will omit this and reinstate the four 
species content when it is appropriate. 

A further coupling of the photon field to the membrane 
coordinates (with coupling constant A) in the form of the 
pull-back of the bulk field strength is also implemented. It 
resembles the Tab term included in Equation ([T]) and rep- 
resents the surface area polarisation energy [I6lll7] . The 
simplest action one can write down for a graphene mem- 
brane interacting with a bulk metric and electromagnetic 
field is then 

Sg^ I d3^(|det[g,b + ^^/^a(A + ^n)^ + AFa,]|)l/^ 

(5) 

where Fab = Fm n daX^^ dbX^ is the pull-back of the bulk 
electromagnetic field strength. We have also allowed for a 
small fermion mass m in the above (the factor of three 
is included for the correct normalisation of the fermion 
mass in the resulting Dirac equation). One should note 
as well that we have two light cone structures; the light 
cone for the electronic properties is given by the Fermi 
velocity, whilst for a bulk vacuum we retain the normal 
Minkowski light cone. As it stands we have not taken into 
account that the membrane is an elastic object. The elastic 
properties can be included by replacing Equation ^ with 



Indeed with a view to including such interactions as in |18l 
I19j , one can write down (for some constant p to be deter- 
mined) 



^iW ^gM^a^rj^M^l^^^ 



(6) 



where the mixed tensor T^ contains the elastic stress ten- 
sor parallel and normal to the membrane and the '^^ (0 
are the fluctuations of the embedding coordinates (the 
phonon modes). In this way one can include the acoustic 
modes (longitudinal and transverse) and their propagation 
in a full manner. 

It is also a simple matter to include a Chern-Simons 
term into the action using the same pull-back principles. 



Sacs = P 



I AaFabe''" 
J Si 



= p f d'^AMF^pdaX^'dbX^d^xPe^"' (7) 
J S3 

which is then added to the dynamical DBI graphene ac- 
tion. In fact as we will see later, we can evaluate this 
coupling constant as a one loop radiatively induced term 
upon integrating out the fermions. 

The coupling constant A, as mentioned previously, has 
the interpretation in the full action as a static surface po- 
larisability due to the effective dipole moments that arise 
from different carbon bonds |16II17] . By performing an ex- 
pansion in Equation ^ to quadratic order in the field 
strengths for static weak fields, one has the schematic re- 
lation A^ ^ otsurface- In coutrast, the vacuum polarisation 
jjMN ^^ (151 results from integrating out the fermions. 
The significance of this coupling is also the way by which 
the photons can interact with the phonon on the graphene 
sheet when considering fluctuating fields. Indeed, from 
Equation ([S]) we see that it is possible to generate three 
and four point vertices that involve both the photon and 
the phonon modes in a self consistent way, similar to what 
is found in [3D] . The main focus for the rest of the paper 
will be on how the XFab term couples to the other fields 
and derive their implications in a few cases of practical 
interest. 



4 Applications and phenomenology of the 
purely bosonic sector: classical aspects 

We can distinguish two cases of practical interest in using 
this action. The first is to perform a derivative expan- 
sion of the action where the membrane degrees of freedom 
are split into background plus fluctuations. The second is 
where one finds a solution of the full field equations for a 
given set of known bulk fields. In both cases we need to be 
able to calculate the determinant of the pull-back fields. 
In this section we will consider only the bosonic mem- 
brane coordinates as the degrees of freedom where we set 
{ip) = 0. Firstly we will consider a toy model of magneto- 
electric media in the bulk to gain some familiarity with 
the pull-back of fields. After this we consider fluctuations 
of the membrane with imposed boundary conditions and 
the effect of the A term. 



4.1 Magneto-electric media in the bulk and a static 
configuration 

An interesting situation to consider is when we have the 
bulk flelds propagating in a background magneto-electric 
media, with an underlying flat metric. For the case when 



James Babington: Graphene membranes and the Dirac-Born-Infeld action 



the permittivity eij{t,x) and permeability fj,ij{t,x) are lo- 
cal functions of spacetime (purely spatial indices i,j — 
1,2,3), we will replace the bulk metric Gmn with both 
electric and magnetic tensors in a simple multiplicative 
way. This naturally produces a local light cone structure 
in the bulk and should properly be thought of as a toy 
model, since we are neglecting all dispersion related prop- 
erties and we putting them in by hand. The pull-back of 
the metric gat in this scheme we choose to be 



gab 



-daX'^dbX" + iefi);.'daX'dbX^. 



(8) 



Consider now a static configuration of a membrane 
immersed in some dielectric with Aj^i — 0. The gauge 
choice of world sheet coordinates are X'~' — ^°, X^ = 
^1, X'^ = ^^, and X^ = Z{^). In addition, we take the 
permeability to be /x^ = 6ij and the permittivity to be 
Cij = diag(l, 1, e(Z)), so that there is only a profile in 
the X^ direction. With this choice the pulled back metric 
becomes 



gab 

-1 + e-^{Z)dtZdtZ diZdtZ d-zZdtZ 

diZdtZ l + t-^{Z)diZdiZ diZd2Z 

d-iZdtZ diZd-iZ 1 + e-\Z)d2Zd2Z 



(9) 



One can now easily calculate the determinant (and square 
root) whereby the action thus becomes 

Sg - Jd'a-^ + e-\Z)[{dtZf - {diZf + [d^Zf] 

+o{z^)Y'^. (10) 

The general scheme of things should then be clear. The 
scalar field will have non-trivial dynamics due to back- 
ground media (in this case dielectric) and the higher deriva- 
tive terms. Since it is a static configuration dtZ = 0. Sup- 
pose it takes the form of an infinite strip of graphene of 
width I that lies in the {X^,X^) plane. The equation of 
motion simplifies when the action is constrained, with the 
area taking some constant value (implying the constancy 
of the effective Lagrangian). Let the model of the per- 
mittivity function profile be given by e{Z) = e ' to be 
definite. Then the equations of motion simplify to 



die 



-Z/2L 



Ci. 



This can be simply integrated to 

Z = -2Lln{ciX^ +C2). 



(11) 



(12) 



What one can draw from this is that the precise shape 
of graphene membranes will be dictated by the nature of 
the static background fields the graphene couples to. As 
alluded to beforehand, this is not a realistic model. Never- 
theless, it is interesting to see that a surface profile can be 
calculated once a background is known. More complicated 
solutions (Catenoids) of this nature can be found in [21 . 



4.2 A graphene drum-skin in an external magnetic field 

The fact that the bulk electromagnetic field couples to the 
membrane in a non-trivial fashion implies that we should 
be able to see its effect in simple phenomenology. To this 
end we calculate the shift of resonant frequencies for a 
graphene sheet that has a disk (drum-skin) configuration. 
This is subjected to a constant magnetic field in a direction 
normal to the graphene disk. 

Consider again the expansion of the action in the static 
gauge detailed subsection 14.11 to quadratic order in the 
fields for only the normal direction (that is the X^ = Z 
direction) fluctuations together with F12 — B^ = B ^ 0. 
We will also include the speed of sound v of the normal 
modes:- 



det(5afc + XFab) = [-1 + v'^Z^ - [diZf - {d2Zf\ 

\2 



+ [-l + i;2z2]A2(Fi2)2. (13) 



This leads to an equation of motion for the transverse 
oscillations 



u^(l + A^B^)Z- V^Z = 0. 



(14) 



One recognizes now a standard drum skin problem. In 
polar coordinates the solution takes the form Z(t, r, 0) = 
e*"*X(r,6') where 



X(r,0) = ^a„J„(A:r)e"^-fc., 

n=0 



(15) 
(16) 



For a drum skin of radius i?, we have the boundary con- 
dition X{R,9) — 0, whereby the frequencies now take on 
a discrete form given by the zeroes of the Bessel function 
Jnikmni) = such that 



w„ 



V2(1 + A2B2)' 



(17) 



In particular, the lowest resonant frequency is shifted to 

Woo « 2.4/(01)2(1 + A2b2). 



5 Fermionic sector applications and 
phenomenology: quantum aspects 

One can see that the fermionic sector is particularly inter- 
esting for the electronic properties of the membranes. We 
have a system in which to study the properties of quantum 
field theories in a low-energy and low-dimension window. 
The DBI action Equation ^ can viewed as a 'machine' 
for producing couplings between the bosonic and fermionic 
sectors, and amongst themselves. As by way of providing 
some examples, we shall investigate the effects of the cou- 
plings between the fermionic terms and the field strength 
XFab in Equation ([5]). In particular, we shall focus atten- 
tion on the scenario where we have static external elec- 
tric and magnetic fields present together with fiuctuating 
fermionic fields. The fermions will be considered in both 
the first and second quantised scenarios. 
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5.1 First quantised fermions in background 
electromagnetic fields: Landau levels 

Before moving on to a specific example let us make a few 
remarks about the mechanics of the DBI action Equa- 
tion ([5|) involving spinors. It is possible to consider the 
fermions on a general curved background (for an overview 
of how curvature affects the electronic properties see [251 
[23] ) and with background electromagnetic fields. The DBI 
action needs to be expanded to quadratic order in the 
fermion fields from which the Dirac equation can be found 
as the leading contribution. To see this, consider the ma- 
trix expansion 

det(M + N) = dct(M) dct(l + M^^N) 

«dct(M)(l + Tr(M-iN)), (18) 

which is true when the absolute values of the matrix ele- 
ments of N are much smaller than those of M. This imme- 
diately yields the massless Dirac equation (for A = 0, m = 
0) on a curved space associated with the pulled back mem- 
brane metric gab in a fully dynamical manner 



1 



= 0. (19) 

Here Fab — [Aj-^c]- It has been used in [Ml[^ to cal- 
culate the energy spectrum of fullerene molecules when 
the curved space is just the two dimensional sphere or a 
spheroid. From now on we consider only flat geometries 
where w^f'^ = 0. 

An interesting application of the model is when there 
are background electromagnetic fields present, arising from 
the A term. The effect of this is to provide a novel type of 
coupling. Heuristically speaking, just as the metric cou- 
ples to the energy momentum tensor of the fluctuating 
fields, the background electromagnetic field will couple to 
the antisymmetric part of the fermion pull-back. 

As an example involving a first quantized field (with 
771 = 0), consider the situation of a plane graphene mem- 
brane (in the x-y plane R^) that is in a constant magnetic 
field B (in the z-direction normal to the graphene plane). 
In this case the expansion is performed about the matrix 



(.ga6 + AFQb)|R,2 




(20) 



To linear order in A we find a modified Dirac equation 

{r''Daii> + \B{r^D2-r2Di)i,) = Q. (21) 

For A = one finds the usual relativistic Landau energy 
levels and wave functions J^ . The energy levels are given 
by E{N) = ±h{y/2eB/c)vFVN, where N = 0,1,2,---, 
whilst the corresponding spinorial eigen-functions are given 
in terms of Hermite polynomials. By performing a first or- 
der perturbative calculation due to the A term, one finds 
a zero energy shift for a given energy level. It is therefore 
necessary to go to at least quadratic order to see its effect 
in connection with the Landau levels. 



5.2 Fermionic fluctuations in background 
electromagnetic fields: induced Chern-Simons terms 
and pair production 

As a further example we will consider the fermionic field 
as a second quantised field. Upon looking at the expansion 
of the DBI action, we see that in addition to the standard 
gauge coupling term eAa4^r°'4^ we have the spin current 
interaction {Xm/6)Fab'ipr'^''ijj. If we compute the one loop 
effective action for these two interactions as in [55], we 
find a three index polarization tensor 7T"'"^(p) given by 



77""'= (p) 






rtr 



T-ia 



(p + kf 



rr" 



TO^ 



m^ 



(22) 



One can extract the low energy behaviour [p — ^ 0) from 
this and indeed it is similar to that of the standard d ~ 
1 + 2 induced Chern-Simons term. The effective action 
reads 



SGCs[A,m,\] 



Xem 

3Tr\Tn\vp 



AaFbcC 



abc 



(23) 



The presence of such terms are of interest when one wishes 
to consider observables such as polarisation rotation ef- 
fects [18 and Casimir-Polder forces [1^. 

As a final example involving the fermion as a quan- 
tum field, one can ask how the Schwinger effect [U] is 
modified due to the presence of the A coupling. Consider 
the situation where we have constant externally applied 
electric field E in the plane of the membrane and the ef- 
fect this has on the fermions when they are taken to be 
second quantized fields. The Schwinger effect for this sys- 
tem has been described in ^^ for a constant electric field, 
together with a proposed experiment to measure the pair 
production in terms of macroscopically observed transient 
currents. The result found there for the rate of production 
per unit area per unit time of fcrmion-hole pairs is 



(eEf/^ 



Y^ 



1 



.3/2, 1/2 A^ „3/2 



exp 



n=l 



eETi 



(24) 



The DBI action Equation ([5]) provides us with a new set 
of couplings. To see the effect of the mass term and the 
polarization term, we must go back to the effective ac- 
tion [5^150] that is obtained from expanding Equation ([S]) 
to quadratic order in the fermionic fields 



lnTy[v4] 



-Tr In 



F'^R, 



F'^{Pa - eAa) - {Xm/6)F^bp^^ _ „, + ^^ 



(25) 



One can evaluate this perturbatively as an expansion in A. 
As for the case of Landau levels, there is no contribution at 
linear order in A. The first contribution arises at quadratic 
order and is found to be 



w 



1+2 



(eE) 



3/2 



2*3/2 1/2 /_^ „3/2 



E 



1 



■ exp 



n7rm2(9-h A2£:2)^;3 



%eEh 



(26) 
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At least for this couphng the effect would be to lessen 
the pair production for increasing electric field strength 
because the polarisation term requires work to be done. 



6 Casimir energies 

One way of evaluating the DBI action is to perform an 
expansion in the fluctuating membrane fields around a 
background. It can be thought of as a 'test particle' or 
probe both in a classical or quantum picture. In particular, 
one may be interested in how the membrane reacts when 
it interacts with other matter. One such scenario is pro- 
vided by Casimir physics [HI] and the subtleties therein. 
For example, in the case of two parallel plates with a gen- 
eral dielectric media in between them, there is a possible 
stress within the dielectric media itself that arises due to 
a particular choice of stress tensor [33]. Thus one would 
expect that a graphene membrane placed in such a media 
as a probe would be subjected to an observable force and 
could be used to distinguish between different theoretical 
proposals. See also [T51IM1IM] for related considerations. 

We will follow here a simpler program and confine our- 
selves to calculating a standard two body Casimir energy 
in the framework of QFT. Specifically, by the considera- 
tion of scattering amplitudes it is possible to calculate a 
static potential between the bodies. 

Suppose we have two infinite parallel planes of graphene 
separated by some distance R, and we would like to know 
what the Casimir energy is between them. It can be found 
by calculating the scattering amplitude for the two - mem- 
brane — >■ two - membrane process. This is what is done 
for the Casimir-Polder potential between two polarisable 
particles based on a phenomenological action [351130] . The 
differences here are that we are dealing with extended (in- 
finite) bodies together with a different form for the effec- 
tive vertices. There are in principal different contributions 
to the static potential; notably the atomic surface polaris- 
ability as well as the contribution due to virtual fermion- 
hole pairs. We shall only consider the contribution that 
arises due to the A term in Equation ([5]). 

The basic scattering amplitude we want to evaluate is 
for two flat and parallel rigid membranes with initial 4- 
momenta kf' and fcf^, and final momenta k[^^ and k2^^ . 
We therefore need to evaluate the four-point correlation 
function (X(1)X(2)X(3)X(4)). The correlation function 
is defined using a two body path integral Zd given by 



Zd = J[dA]exp Uj F^Ia\ Wd[Al 

2 / 2 

y^d[A] = [tl [dX^di^'d^'] exp j J2 iS, 



(27) 



(/' . (28) 



In this setup, all of the fields are taken to be full quantum 
fields. We are therefore interested in obtaining another 
effective action and thereby vertex, where we know how 
the field strength Fmn couples to the membrane X*^. 
The DBI action Equation ^ already contains such a set 
of terms given by the A term. If we expand out the bosonic 



part to second order in A, together with the gauge choice 
of coordinates Equation ^ (taking the elastic matrix to 
be Tj^i — 6^j for simplicity) one can obtain the necessary 
couplings 



S[X,F]^ Jd^^[\dctigab + XFa, 



11/2 



1 - \da<l'^' da'^'"' - ^Fat^F'^" 



F„ 



Fmn {5 a 



da<P''){5'^ +d,<P'') 



(29) 
(30) 



In this way one generates both the standard projection 
of the the full field strength onto the membrane, a cubic 
coupling and most importantly a quartic coupling of the 
form 



S^tl-Jd'(X'F'{d<Pf 



(31) 



With this effective vertex it is a simple matter to cal- 
culate the interaction energy between the two graphene 
membranes. The four-point function becomes 

{X{1)X{2)X{3)X{4)) ^ X" f d^^id^UX{l)dX) 



x{X{2)dX){Xi3)dX){X{4)dX) 
x{F{X{^,))F{X{^,))) 
x(f(X(6))^(X(ei))), (32) 

where ^i and ^2 are the world-sheet coordinates on mem- 
brane one and two respectively. From the known form of 
the time ordered photon propagator, the correlation func- 
tion (FF) is given by 



{F{X{^,))F{X{^,))) 



1 



[R' + \ci~^2\r 



(33) 



The remaining two point function of the scalars just de- 
scribe the in and out legs. It is possible now to calculate 
the corresponding S-matrix and the static potential. By 
doing the time integration and then the integration over 
all the two dimensional membrane space, one derives the 
static potential density. This can clearly be seen to behave 
like V ^ X'^/R^. It is interesting to see a different power 
law behaviour that arises due to such a phenomenological 
interaction term. 



7 Conclusions 

The intent of this paper has been to set up a simple phe- 
nomenological scheme to address the dynamics of graphene 
sheets as a whole. This includes both the overall dynamical 
structure of the graphene lattice approximated as a con- 
tinuum and their electronic properties as described by the 
inclusion of world sheet fermions. These can be naturally 
coupled to the bulk electromagnetic field and different 
schemes exist with which to evaluate the resulting action. 
One notable feature is the way in which the fermion spino- 
rial fields live on the graphene membrane and couple to 
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bulk fields which in principle provide a consistent scheme 
for calculating bulk-membrane quantities, e.g. graphene 
sheet reflection coefficients in a dielectric. The inclusion 
of a pull-backed electromagnetic field into the action (the 
A term) as a way of including the static polarisation of 
the membrane has consequences when we consider field 
fluctuations. Upon expansion we find cross couplings that 
produce additional perturbative terms, that in principle 
have observable consequences. 

An important and obvious next step is to consider the 
effect of real materials in the bulk and finite temperature. 
All the fluctuating physical fields need to be Fourier trans- 
formed in time as well as the bulk background fields. It 
would also be interesting to try and apply the DBI formal- 
ism directly to scattering experiments in the low energy 
sector. The calculation of photon scattering amplitudes 
from a graphene sheet should give the same information as 
the reflection coefficients obtained from matching condi- 
tions. Pursuing a S-matrix programme would be a worth- 
while effort to better understand the associated scattering 
phenomenology. 
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